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Abstract 

The notion of s-points has been introduced by the authors (SIAM 
JMA, 39 (2008), 1821-1850) in connection with the control prob- 
lem for the dynamical system governed by the 3d acoustical equation 
uu — An + qu = with a real potential q G Cq°(K 3 ) and controlled by 
incoming spherical waves. In the generic case, this system is control- 
lable in the relevant sense, whereas a G M 3 is called a s -point (we write 
a G T q ) if the system with the shifted potential q a = q(- — a) is not 
controllable. Such a lack of controllability is related to the subtle phys- 
ical effect: in the system with the potential q a there exist the finite 
energy waves vanishing in the past and future cones simultaneously. 
The subject of the paper is the set T q : we reveal its relation to the 
factorization of the S'-matrix, connections with the discrete spectrum 
of the Schrodinger operator — A + q and the jet degeneration of the 
polynomially growing solutions to the equation ( — A + q)p = 0. 



Introduction 

0.1 Dynamical system 

An acoustical scattering problem is the system of the form 

u tt - Au + qu = 0, (x,t) G M 3 x (-00,00) (0.1) 
u| w< _ t =0, t<0 (0.2) 
lim su((s + t)u),-s) = f(r,uj), (r, ui) G [0, 00) x S 2 , (0.3) 

where u = u^(x,t) is a solution (wave), q = q(x) is a real valued smooth^] 
compactly supported function (potential), and / G T := L 2 ([0, 00); L 2 (S 2 )) 



1 everywhere in the parer "smooth" means C°°-smooth. 
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is a control. With the system one associates the control operator W : T — >■ 
H := L 2 (IR 3 ), M 7 / := ( • , 0) and the subspaces: U := RanVF {reachable 
set), T> := T-LQlA [defect subspace), M := KerVF C J 7 [null control subspace). 
The relations < dimT> = dim M < oo hold [2]. If T> — {0}, the system 
flOIIJ)-(IE9) is said to be controllable H; the case of T> ^ {0} is referred to as 
a lack of controllability. 

0. 2 s-points 

As is shown in [2], the case dim 22 ^ is realizable and corresponds to the 
curious effect: for / G M, the function w^(x,t) := f Q (x , s) ds is a finite 
energy solution of ( 10. ip . which satisfies w'{ - , —t) — w*(- ,t) and 

W / | larKltl = , (0.4) 

1. e., describes the wave vanishing in the past and future cones simultaneously. 
This wave comes from infinity, stops at the moment t — 0, and then return 
back to infinity along the same trajectory. Such a behavior motivates to 
call it a reversing wave and regard the coordinate system origin x = as a 
point, which is able to stop incoming waves ('stop point'). Looking for such 
points in the space, we change the coordinates x > x — a and deal with 
the system f lU.ip ( IDTSj) with the shifted potential q a := q(x — a), the objects 
corresponding to q a being labelled with the subscript a. An a G M 3 is said 
to be a s-point of the potential q if T> a ^ {0} holds. By T q we denote the set 
of such points. 

To the best of our knowledge, s-points is something new in the acoustical 
scattering and the set T q is worth studying for its own sake. The goal of our 
paper is to originate such a study: we reveal certain relations between T q 
and the known objects of the scattering theory. 

0.3 Results 

Evolution of the system (I0.ip - fl0.3p is governed by the Schrodinger operator 
H := —A + q in TL, DomiJ = H 2 (M. 3 ). This operator may have a finite 
discrete spectrum adisc{H) C (— oo, 0] (see, e.g., [8]) and in the paper, for the 
sake of simplicity, we accept 

2 Note that the unperturbed system (with q = 0) is controllable. 
3 We mean w{( ■ , 0) = 0. 
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Convention 1 The potential q is such that the equation 

(-A + q) tp = in M 3 (0.5) 

has not nonzero solutions, which satisfy if(x) — > as \x\ — > 0. 

In physical terms, this means that the hamiltonian H has neither bound nor 
semi-bound state at the zero energy level (see, e.g., [6]). If exists, such a level 
can be removed by arbitrarily small perturbation of the potential q, so that 
the convention is not restrictive and we deal with a generic case. 
Our results are the following. 

• Let S a {k), k G IR be the Schrodinger S- matrix of the potential q a . In 
the series of papers cited in his book [6], R.Newton introduced the 
Riemann-Hilbert factorization of the form 

S a (k) = U-(k)[S-(k)S+(k)]Ut(k), keR, (0.6) 

where IT^ are the Blaschke-type rational operator-valued functions, 
are the operator-valued functions holomorphic and boundedly invert- 
ible in the half-planes {k G C | ± Qk > 0} respectively. This represen- 
tation is used for solving the inverse problem that is to determine the 
potential q from the scattering data: it is shown that the knowledge of 
the family {^jagRS enables one to recover q. In the mean time, in the 
mentioned papers it is assumed that such a factorization is realizable 
for any a G M 3 . We show that for a G T q the representation W.6\) fails, 
i.e., the factorization is impossible^. 

• In [2] we suggested that the presence of s-points is connected with the 
discrete spectrum of the operator H: the conjecture was that T q ^ 
is equivalent to a^ lsc {H) ^ 0. Here this conjecture is partially justified 
as follows. 

Let a point o 6 I 3 , an integer m > 1, and a function h G % be such 
that 

(-A + q) m h = inR 3 \{a} (0.7) 

holds, whereas (—A + q) m ~ h does not vanish identically. As is shown 
in [2], such an a is necessarily an s-point; by T q m C T q we denote 
the set of these points and specify m as the order of a. We prove that 
Tq 1 7^ is equivalent to a^ sc {H) ^ 0. 

4 However, we do not claim that this fact cancels the determination of the potential by 
R.Newton's procedure. 
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• We show that the points a 6 T g fin := U m >iT 9 m can be specified 
as the jet degeneration points of polynomially growing solutions to the 
equations (—A + q) p = in R 3 . By this, typically the set T g fin consists 
of s-surfaces. If the potential is radially symmetric, i.e., q = q(\x\), the 
s-surfaces are spheres and their position in the space can be studied in 
more detail. 

The work is supported by the RFBR grants 08-01-00511 and NSh-4210.2010.1. 

1 s-points and R.Newton's factorization 

1.1 Scattering matrix 

Recall the basic definitions and facts (see, e.g., [6]). 

The Schrodinger scattering operator of the pair Hq — — A, H = —A + q 
is S : U -+ H, S := W^W_, where W± := s - lim Moa e' itH e itH ° are the 
wave operators. 

3 

Let (Fy) (p) = (27r)~ 5 J r3 e ipx y(x) dx be the Fourier transform. The ope- 
rator S := FSF" 1 is of the form 

(&) (p) = (S(k) [v(k, ■ )]) (0) , 

where v(fc,o;) := v(ku), k := |p|, 6> := |y, uj,9 £ S 2 := {p G 1R 3 | |p| = 1}, 
whereas S(k) is an operator acting in L2(S 2 ) on the angle variable by the 
rule 

(S(k)g)(6):= [ a(0,a>',k)g(u)du, 
Js 2 

depending on /c as a parameter and being called the S-matrix of the potential 
q. So, we regard the S-matrix as an operator-valued function of k > 0, the 
values being taken in the bounded operator algebra B(L 2 (S 2 )). The kernel 
s(9, to; k) is a distribution of the form 8(9 — u) +s"(0, oo; k) with a smooth so 
that each S(k) is of the form "identity + compact operator". Moreover, each 
S(k) is a unitary operator. Also, the well-known high energy asymptotic 
S(k) — > I as k — > oo holds. 

5 Presumably, T q n exhausts T g but it is still a conjecture. 
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In what follows, we assume S(k) extended to k < by S(k) := S*(—k). 
The index of the S-matrix is defined by 

ind5:=i[lBdfit5(Jb)]|^; 

by the Levinson theorem it is equal to the total multiplicity of the discrete 
spectrum a<n sc (H). The representation 

S(k) = U-{k)S red {k)U + {k) , keR (1.1) 

is valid, where S Ted is the so-called reduced S-matrix, which satisfies ind>S red = 
0, and II is a B (L2(>S 2 ))-valued function of the form 

n± (*» = n( i+ i3^^) • k£m 

with the certain projections = ( • , ijj^) L r S 2\ ip^i w h ere — "0~(~ w )) 

and k n > such that — k 2 n 6 cr disc (i?) (see [6]). 



1.2 Factorization by R.Newton 

In solving the scattering inverse problem by R.Newton, the following Riemann- 
Hilbert type representation of the reduced S-matrix plays the key role: 

S red (k) = S-(k)S + (k) , feel, (1.2) 

where S are the B (/^(S^-valued operator-functions, which are holomor- 
phic, bounded, and boundedly invertible in the complex half-planes C ± = 
{k G C | ± Qk > 0} respectively. In more detail, to recover the potential q 
via its S-matrix, one needs 

• for a fixed a 6 K 3 , to find the S-matrix S a of the shifted potential q a 
that can be done in a simple and explicit way 

• to determine S red a from (11. ip and represent 5' red a = S~S+ by (II. 2p 

• by varying a, to collect the families {^} oeR3 ; each of them determines 
the potential q. 
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However, in [6] the author assumes that the second step can be fulfilled for 
all a G M 3 . The following result shows that such an assumption may be 
invalid. 

Theorem 1 If a G T q then the representation M. 6 ^) for S Tcd a does not hold, 
i.e., the factorization S Ted a = S~S^ is impossible. 

We postpone the proof till section 1.4 and begin with preliminaries concerning 
to the well-known facts of the Lax-Phillips theory. 

1.3 On Lax-Phillips scheme 

The Cauchy problem for the acoustical equation is the system 

v tt - Av + qv = 0, (x,t) G M 3 x (-00,00) (1.3) 
v \ t=0 = <p, v t U=o= V> in (1.4) 

with the finite energy data d := {tp, ip} G D := // 1 (M 3 ) x L 2 (M 3 ); the solution 
is denoted by v = v d (x,t). A peculiarity of the acoustical scattering is that 
the energy form 

E[d,d'} := / + V<f/ + q<p<f/, 

which the set of data is equipped with, is indefinite. However, the form 
turns out to be positive definite on the absolutely continuous subspace D ac : = 
DQ E D disc , where D disc is the finite dimensional subspace spanned on the data 
{(fk, ±k(pk} such that — k 2 G aai SC {H) and Hipk = —k 2 tpk- By P ac we denote 
the projection on the first summand of the decomposition D = D ac ®e -Ddisc- 
In the framework of the Lax-Phillips theory [5] , there is a certain freedom 
in the choice of the pair of the incoming and outgoing subspaces D . Once 
such a choice made, for each d G D the corresponding incoming and outgoing 
spectral representatives d T , which are the 5* 2 - valued functions of k G K of the 
class L 2 (M; L 2 (S 2 )), do appear. In the spectral representation, the scattering 
process is described by the Lax-Phillips scattering operator S : gL i-> d + , 
which acts in L 2 (R; L 2 (S 2 )) by 

(Sg)(k) = S(k)g(k), fcGK, 

where S(-) is a B (L 2 (5' 2 ))-valued operator-function (the Lax-Phillips S- 
matrix). Two possible variants of the choice are the following. 
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1. Asskr n a d G D to the subspace Dq C -D if v \\x\<—t — for t < and 
the subspace C D if v d || x |<t= for t > 0. The subspaces 

D T := P &C DT 

constitute an incoming/outgoing pair [5]. By S and S(-) we denote 
the corresponding scattering operator and its S'-matrix. 

2. One can reduce the incoming/outgoing subspaces to the smaller ones 

Dl A := Djn^C^. 

By S red and S red (-) we denote the corresponding operator and S'- 
matrix. 

The important fact, which relates the quantum scattering and acoustical 
scattering objects and follows from the above accepted definitions, is that 
the equalities 

S(k) = S{k) , S Icd (k) = S vcd (k) , keR (1.5) 

hold (see [5], chapter VI, part 2, the relation (3.4)). Note in addition that 
( TT3|) and drTJ imply 

S(k) = U-(k)S red {k)U + {k) , keR, 

whereas the factors n T ( ■ ) can be also interpreted in terms of the Lax- 
Phillips theory as S-matrices corresponding to the certain choice of the in- 
coming/outgoing subspaces. 
Let 

L 2 (R; L 2 {S 2 )) =U~ ®U + (1.6) 

be the decomposition on the Hardy subspaces which consist of functions 
holomorphic in {k G C | =F Qk > 0} respectively; by P T we denote the 
projections on H T . As is known, if d belongs to the incoming (outgoing) 
subspace then its spectral representative lies in %~ {H + ). Representing the 
scattering operator in the matrix form in accordance with (jl.6p . one has 

/ p— cred p— p— cred pH 
r>red _ I r »J •* X. I (l 71 

1 "P+^rcdjj— p+gredp+ I > ' 
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1.4 Proof of Theorem [T] 

Assume that a G IR 3 is a s-point of the potential q. As is evident, by shifting 
the coordinate system, one can provide a = 0. So, let G T r 

Let wf be a reversing wave (see section 0.2), d := {u>^( ■ , 0), 0} its Cauchy 
data. By ( 10. 4(1 . one has <io G -D - fl D + ; moreover, the results [2j on the 
stability of trajectories easily imply do G D~ ed PI -D+ d . The latter leads to 
(io G H~ and c£q G H + . Therefore, 

S red dQ = g?q~ G 

and, applying the projection P~ , one gets P~S rcd dQ = P~ S rcd P~d = 0. 
Thus, we see that Ker P~S red P- ^ {0}, i.e., the block P~S rcd P- in (TUD is 
not invertible. 

By (ll.5p . the block p-s ved P~ is also not invertible. However, as is well 
known, such an invertibility is the necessary condition that provides the 
Riemann-Hilbert factorization (11.21) (see, e.g., [6]). Hence, in the case of 
a G T g , to represent S red a in the form (II. 2p is impossible. □ 

In addition, note the following. Although the representation (II. 2p fails 
for a 6 T 9 , by very general results of the Riemann-Hilbert theory one can 
factorize the reduced S-matrix as 

S^ d a (k) = fl-(k) [S-(k)S+(k)] fl+(k) , keR (1.8) 

with the Blaschke-type operator functions Il a , A*, which have the certain 
complex poles, and holomorphic and boundedly invertible in {k G C | ± 
Qk > 0} [3j. Moreover, it is the presence of the Blaschke factors, which 
renders the representation (11.21) impossible. In this connection, the intriguing 
questions arises: What is the physical meaning of the factorization (II. 8p and 
the corresponding Blaschke poles? Can one characterize them in terms of 
the Lax-Phillips theory? The questions are open. 

2 s-points and discrete spectrum 
2.1 Theorem 2 

Assume that for a point a G M 3 there exist a function h G % and an integer 
k > 1 such that 

(-A + q) k h = inM 3 \{a} 
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and let m be the minimal of such fc's. As is shown in [2J, this a is necessarily 
an s-point, whereas m is specified as the order of a. By T q m C T 9 we denote 
the set of s-points of the order m. 

Remark In the light of these definitions, Convention [1] is motivated as 
follows. If G (Jdisc(H) and (p G H is a zero energy level eigenfunc- 
tion, then (—A + q) fo = everywhere in R 3 and we are forced to accept 
Tg 1 = T q = R 3 . The convention excludes such a degeneration. 

The following result partially explains why s-points do appear: one of the 
reasons is the presence of the discrete spectrum of the operator H. 

Theorem 2 T^ 1 ^ holds if and only if o~disc(H) ^ 0. 

The rest of section 2 is the proof of the theorem. 



2.2 Green function 

We define the Green function of the operator H ClS db Li, oc (R 3 ) -solution of the 
integral equation 

1 f l( S ) n!„ .A j„ _ ^ Tn.3 



G(x,2/) = — = 1 - / ; w a; G R , (2.1) 

where y G R 3 is a parameter. This definition is correct since this equation 
is uniquely solvable. Indeed, otherwise there is a nonzero solution ip of the 
homogeneous equation 

<p(x)=-[ t , , q(s) ds , xGR 3 . 
J RS An\x - s\ 

As is easy to see, such a solution satisfies ( 10. 5ft and vanishes at infinity that 
is forbidden by Convention [U To establish the solvability, one can reduce 
the equation to a large enough ball B D supp q, prove the solvability there, 
and then extend the solution to R 3 \I? h.s. of (TJU). 

Also, the integral equation implies the symmetry property G(x, y) = 
G(y,x), x 7^ y. From ( 12. lft one derives that in the sense of distributions 
on C^°(R 3 ) the function G satisfies 

(-A + q)G(-,y) = 8 v , (2.2) 

where 5 y is the Dirac measure supported in y. 
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Lemma 1 The asymptotic 

G(x,y) = Pr\+o(r^) , yeR 3 (2.3) 
\x\-hx> An\x\ \ \x\ 2 J 

holds, where $ is a smooth function obeying 

(-A + q) $ = in R 3 (2.4) 

<&(x) ->• 1 as |x| -)■ oo , (2.5) 

the asymptotics \2.3\) and A2.5\) being uniform w.r.t. A G S* 2 . 

Proof Fix x in the equation (12. ip and tend \y\ — > oo. Using the symmetry 
property and looking for the solution in the form (12. 3p . we have 



m +0 (A.) = _^ +0 (-L) _ f gW 

that easily leads to the integral equation 



47r|y| \ j/ 



ds 



$( x ) = l-/ $( a ) rfg x g R 3 . (2.6) 

J M3 4vr|x - s\ 

The latter equation implies ( 12.41) . ( 12.5j) . and justifies ( I2.3p . □ 
Now we are ready for proving Theorem [2j 



2.3 Necessity 

Let 0"di SC (-fO 7^ and = — k%x f° r = m f cr disc(^), so that % is 
a ground state of the operator H. As is well known, — & 2 is an ordinary 
eigenvalue of H, whereas the eigenfunction x behaves as x( x ) ~ e - fc oM 

and can be chosen positive: x > everywhere in R 3 . Integrating by parts, 
we have 

= (see (|23J|) = / x(*0 [(-A + g) $] (z) dx = -fcj} / ds. 

Hence, $ has to change sign and there is an a G R 3 such that $(a) = 0. As 
result, we conclude that G( ■ , a) G % (see (12.31) ). whereas (—A + q) G( ■ , a) = 
in R 3 \ {a} (see Thus, a G T ? a and, hence, T^ 1 ^ 0. 
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2.4 Sufficiency 

The sufficiency will be proved by contradiction. Assume that o"di SC (-^) = 
but T^ 1 7^ and, hence, there is a point oGl 3 and a function ft E H such 
that (-A + q) ft = holds in M 3 \ {a}. 

Show that $(a) = 0. Indeed, considering (—A + g) ft as a distribution on 
C£° (M 3 ), we see that it is supported at x = a. Such a distribution has to be 
a linear combination of the Dirac measure derivatives: 

N N 

(-A + q)h=J2 a i D * S " = E n > { l)Jl) '^ < 

b'l=o I i |=o 

where j = {jx,j2,h} is a multi- index, |j| = j 1 + j 2 + j 3 and D 3 X = d^d^d 3 ^ 
is a differentiation, cu^ are constants (see, e.g, j9]). Therefore, by fl 2 . 2 j) . the 
function /i has to be of the form 

N 

h(x) = a J {-l) ljl DiG{x,a) + h{x) , (2.7) 
b'l=o 

with a h satisfying — (A + q)h = as a distribution. Hence, ft, is a smooth 
function and 

- (A + q)h = in M 3 (2.8) 

holds in the classical sense. As can be shown from (12. ip . for |j| > 1 the 
derivatives D 3 a G( ■ , a) are not square-summable near x = a. Therefore, the 
only option for the sum in (12.71) (and, hence, for h) to belong to H (to be 
square-summable near x = a) is otj = for \j\ > 1 and we get 

/i(sc) = a C(2;, a) + ft(x) in M 3 \{0}. 

For large enough \x\ > diam supp q, we have Ah = and know that h E~H. 
Such a function vanishes as |a;| — > oo. Hence, by (12.31) . the function h also 
tends to zero as \x\ — > oo and, in the same time, satisfies (12.81) . By Convention 
1, one has h = 0. Thus, ft is proportional to the Green function G( ■ ,a), 
whereas G(- ,a) eU implies $(o) = by (J2T3]1 . 

Now we apply the perturbation theory arguments. Consider an operator 
family {—A + £(?} £e [o,i]; let $ £ be the corresponding analog of the function 
$ = $ x that appears in the asymptotic (I2.3p . Note that for e = one has 
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$° = 1. As is easy to see from the integral equation ( 12. ip . for small s the 
inequality $ £ ( ■ ) > holds everywhere in R 3 , whereas 



<& £ (x) — >• 1 as \x\ — > oo 



(2.9) 



uniformly w.r.t. e. 

Since $ £ depends on e continuously 0, there is the minimal sq G (0, 1] 
such that the function $ eo satisfies <3> £0 > and, in the mean time, has zeros. 
Shifting (if necessary) the origin of coordinates, assume that $ £0 (0) = and 
recall that $ £0 satisfies (-A + e q) $ £0 = in 1R 3 . 

Fix a positive r; let B r := {x G R 3 | |x| < r} and S r := dB r . Integration 
by parts implies 

9(?) * e " K) (i - 9 d ^ 1 A *" K) (ili - 



From the other hand, we have 

r 

< max \ q\r I U(r) dr < 

Br J 



< max \q\ j dr r 2 

Br 



1 1 

r r 



<$> £0 (t9) d6 



Si 



max |g| 

Br 



max [/ 

[0,r] 



and arrive at the estimate 

U(r) < cr 2 max U(-). 

[0,r] 

Choosing r small enough such that £7(r) = max[/( ■ ) we see that this esti- 

[0,r] 

mate is possible only if U{r\) vanishes identically for ri < r. Since $ eo > 0, 
this means that $ eo = in a small ball B r . 

So, $ £0 is a solution of an elliptic equation vanishing in a ball. By the 
Landis uniqueness theorem [5], such a solution vanishes everywhere and we 
have $ £ ° = in R 3 that contradicts to (12T91 . □ 



6 in fact, analytically 
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3 s-points and degeneration of jets 
3.1 Jets and polynomials 

Recall what a jet is. For a fixed a G R 3 and an integer k > 0, we say the 
smooth functions u and v to be equivalent (and write u ~ v) if — = 
o (|x — a\ k ) as x — > a. With respect to the relation ~, the set C°°(R 3 ) decays 
on the equivalence classes, whereas the class j k [u) := {v \ v ~ u} is called the 
jet of u (at the point a, of the order k). Setting up ajaftt] +/3j%[u] := j k [au + 
f3v] for a, (5 G R, one makes the set of jets into a linear space J k [C°°(1R 3 )]. 

As is easy to recognize, the jet can be identified with the collection 
of the derivatives {D J x u(a)}y^ =0 , whereas dimjj [C°°(1R 3 )] is equal to the 
number d k of such derivatives: d° = 1, d 1 = 1 + 3 = 4, d 2 = 1 + 3 + 6 = 
10, ... ,d k = i(Jfc + l)(jfe + 2)(fc + 3). 

Let C C C°°(1R 3 ) be a linear set; introduce a subspace 

== W I « G C] C J* [C°°(R 3 )] 

and denote by its dimension. 

Fix an integer / > 0; the functions belonging to a linear set 

V\ := {p G C°°(M 3 ) | (-A + g)p = inR 3 , |p(x)| < const (1 + |a;|)'} . 

are said to be q-harmonic polynomials of the order < I. The existence of such 
functions is established in a standard way: one puts p = pi + w, where pi is a 
harmonic polynomial (see, e.g., [9]), derives the relevant integral equation for 
w of the form analogous to (12.11) . and proves its solvability and uniqueness of 
the solution by the same arguments as for the equation (12.11) . i.e., referring 
to Convention [TJ For instance, the function $ determined by (12.41) . (12.51) is 
an element of V®. Also, as is easy to verify, V l q C V l q for I < I', and the 
relations 

dimPj = dimV l = (I + l) 2 , d k a [V k q }<d k a {V k } = {k + lf (3.1) 

hold. 

The main result of this section is 

Theorem 3 For m = 2,3, ... , the inclusion a G T q m is equivalent to the 
relation 

dl m - 2 [P^ m - 2 ] < (2m- l) 2 . (3.2) 
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Since in the generic case one has d 2m ~ 2 [P 2m ~ 2 ] = (2m — l) 2 (see (13. ip for 
k = 2m — 2), it is reasonable to refer to ( 13. 2 p as a jet degeneration of p 2m ~ 2 
at the point x — a. Also, the Theorem [2] can be interpreted in the same 
terms: $(a) = means that dim jo fPj] = < 1. 

The proof is postponed till section 3.3. For the sake of simplicity, it will 
be demonstrated for the case m = 2; the way to treat the general case will 
be clear. 



3.2 Green function 

Here we deal with the g-biharmonic equation (— A + q) 2 u = h. For x, y G M 3 , 
(x, y) := Xiyi + x 2 y2 + x 3 y 3 is the inner product. 
Unperturbed case For q = 0, we set 

d t \ \ x ~y\ 

B {x,y) := 

and, in view of the well-known relation 

(-A) 2 B (-,y) = 6 y , (3.3) 

refer to Bo as the unperturbed Green function. The following is a standard 
way to derive the asymptotic of Bo as \x\ — > oo: 

x-y\ = (\x\ 2 -2(x,y) + \y\ 2 ) h = 



x 



x 



X 



X 



X 



, >l 2 2 / x 



2 \\x\ 2 \x\ \\x\' " / J 8 \\x\ 2 l x l \ I 2 - 
M 2 1 / x _.\ 1 / x „\ 2 , ,-wl 



1 + W"H\R7"W\R7 +0(W ' 
+ J^! _/* „\ _J_ /«, „V + 0(W -») = 



2|x| \ \x\ I 2\x\ \ \x 



1 2/ 1 2 1 \ - 1 x - 1 x - 2 2 

II; j^j -I 
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where i and j run over 1,2,3. To clarify the structure of this expression, one 
can use the identity 



E2 2 
Xi Vi 



4) {yl - yl) + gfa? + 4 - (vl + vl ~ 2 vl) + 3M 2 \v\ 



and write the result in the form 

12 



\x - y\ 



\x\ + 



\v\ 



3 s 



\x\ 



\x\ 



ViVj 



l 

4 
+ 



>i 2 „2 

1 2 



(I/? - vl) ~ t; 
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,-y>2 _|_ ™ 2 



2xi 



x 



(y 2 1 +y!-2yi) 



x\ 



(3.4) 



As is seen now, this is an expansion over the spherical harmonics: the terms 
into the first curly braces do not depend on angle variables, the ones in 
the second and third braces are proportional to ^/"(rfr) and ^"(jfj) respec- 
tive^. One more peculiarity of such a representation is that among the 
^/-dependent coefficients 



$o = l, q°=\y\ 2 , 
$° = Vl , $° = y 2 , & A 



1/3 , 



$° = y m , $° 6 = y 2 y 3 , $° = vm , $° = vl - vl > $9 = vl + vl - 2 vl 

(3.5) 

there are nine harmonic functions $5) • • • > ^9 an d one biharmonic 

Remark 1 So, the leading part of the unperturbed Green function asymptotic 
\3.$ is of the following structure: 

• the term harmonic w.r.t. x has the coefficient \1/ biharmonic w.r.t. y 

• the other terms are biharmonic w.r.t. x and have the coefficients $° 
harmonic w.r.t. y. 



7 the harmonics arc selected by the square brackets 
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Quite analogous picture will be observed in the perturbed case. 

Perturbed case For q ^ 0, we define the (perturbed) Green function as a 
L^ C (R 3 ) -solution of the equation 

B(x,y) = - \ x ~ y \ + f LB{t,y)dt, xeM 3 , (3.6) 



57T J R3 87T 

where i/ 6 R 3 is a parameter and L := (—A + g) 2 — (—A) 2 is a second-order 
differential operator^ such that supp Lu C supp q fl supp if. The following 
result shows that this definition is correct. 

Lemma 2 The equation A3. 6]) is uniquely solvable. 

Proof Let the homogeneous equation 

= [ ^ -mt)dt 

have a nonzero solution if). Such a solution satisfies 

(-A + qfip = in R 3 ; 
by (13.41) . for large \x\ it behaves as 

3 

ip(x) = a\x\ + y PiT-r + Tpix) , 

|x|-Kxj \X\ 

8=1 1 1 

where if) = — > 0. 

Assume that if) is g-harmonic: (—A + q)if) = in M 3 . If so, we have 
Aif) = for large |x| that implies a = = and if) = if). Hence, if) is a 
g-harmonic function vanishing as |x| — > 00 that is forbidden by Convention 

m 

Thus, ip := (—A + q) if) is a nonzero function, which satisfies (—A + q) tp = 
in M 3 and 



<p(x) 



- a A\x\ - J2 A A fr + A ^ = ( n) ^ 

i=l "' VFl/ 



3 in the integral, L acts on the variable t 
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that is impossible in view of Convention [TJ So, assuming ip ^ we arrive at 
a contradiction. □ 



As easily follows from ( 13. 6p . the relation 

(-A + q) 2 B(.,y) = 8 y 



(3.7) 



analogous to (13. 3p is valid. Also, the symmetry property B(x,y) = B(y,x) 
holds. 

The way to derive the asymptotic of B analogous to (13 .4p is also quite 
standard. Namely, we represent 



B(x,y) =-— A(x,y) + 0(\x\ 



57T 



(3.8) 



with an ansatz of the forrr|§ 

A(x,y) = \ \x\Mv) + 



*(y) 

3\x\ 



\x\ 



1 

\x\ 

1 

4 



My) + 

X\X2 



x 2 



X 





1 2 



x\ 

My) + 



My) + 



£3 

\x\ 



My) - 



x 2 x 3 


My) + 


X1X3 






. \ x \ 2 . 




\x 2 



My) + 



\x\ 



My) - T2 



/y>2 I ,-yi 2 Q™2 

Ju J I tC' 2 .X < g 



$9(2/) 



(3.9) 



tend |y| — )■ 00 in (13. 6p . using the symmetry property substitute this repre- 
sentation to ( 13. 6p . and by equaling the coefficients at the proper spherical 
harmonics, eventually get the 'transport equations' of the form 



= $°(x) + 



I a; — 1\ 

8tt 
|x — t| 

R3 87T 



L$i(t)dt, 
LV(t)dt. 



% = h 



9 



(3.10) 
(3.11) 



By the same arguments, which provide the unique solvability of (13. 6p . these 
equations are also solvable uniquely. 

As is easy to see from (I3.10P and (13.111) . the solutions ^ and $j are 
q-biharmonic. Moreover, 



5 compare with (13. 4j) ! 
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Lemma 3 The functions $1, . . . , $9 are q-harmonic polynomials of the or- 
der < 2 and constitute a basis in Vg. 

Proof 

Step 1: q-harmonicity Take a / G C^°(IR 3 ), denote H : = (—A + q) and 
put 

F(x):= [ B(x,y)Hf(y)dy, x G M 3 . (3.12) 
Applying and taking into account ( 13. 7p . we get H 2 F = Hf and, hence, 

H (HF-f) =0 in M 3 . 

As is easily seen from (I3.12p and (13. 8p . (13. 9p . the function HF vanishes as 
\x\ 00. Hence, the same is valid for HF — f. In the mean time, the 
latter difference is annulated by H and, whence, it has to vanish identically 
by Convention [TJ Thus, we arrive at HF — f = 0, i e., 

(-A + q(x)) f B(x,y)Hf(y)dy = f(x) , x G M 3 . (3.13) 

Now, tending \x\ — > 00 and taking into account the structure of the 
asymptotic (13. 8p . (13. 9p . it is easy to see that the asymptotic of the l.h.s. 
contains certain linearly independent harmonic (w.r.t. x) termJ"1 with the 
coefficients of the form J R3 $«(?/) Hf(y) dy. In the mean time, the r.h.s. in 
the equality (13 . 1 3[) is compactly supported and, hence, the equality forces all 
harmonic terms to vanish identically. This follows to 

0= / <S> i (y)Hf{y)dy= [ H^(y)f(y)dy 

that is equivalent to = in I 3 by arbitrariness of /. So, $j are q- 
harmonic. 

Step 2: integral equations The functions $j satisfy the integral equations, 
which are different from and more informative than (I3.10p . These equations 
can be derived as follows. 

So, H is understood as a differential expression that can be applied to any smooth 
enough function not necessarily belonging to %. 

n for large enough \x\, 'harmonic' and 'q-harmonic' is the same 
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By (E2D, for any / G G£°(R 3 ) one has the relation 

(-A + q(x)) [ G(x, y) f(y) dy = f(x) , x G R 3 . 
Jm. 3 

Comparing it with (13. 13|) we get 

G(x,y)f(y)dy = [ B(x,y) Hf(y)dy , 
Jr 3 

whereas integration by parts easily implies 

(-A + q(y))B(x,y) = G{x,y) , x,y Gl 3 , x^y. (3.14) 
In the unperturbed case, for the Green functions B and G we have just 



-A 



f ~ y\ 

8tt 



4ir\x — y\ 



x,yeR 3 , x^ y. (3.15) 



These relations enable one to get the asymptotics for G and Go through 
(13.81) and (13. 4p respectively, the coefficients of the asymptotics being expressed 
through \l/ and ^° . Surely, the asymptotic of G obtained in this way is 
just a more detailed version of H 2 . 3 f) . which takes into account the structure 
of the lower order terms. Also, the comparison of the coefficients at A in 
the detailed asymptotic and (12.31) easily implies the important relation 

<%) = -±(A + q(y))y(y), (3.16) 
o 

which will be required later on. 

Thereafter, the trick, which was used for derivation of the equation ( 12. 6p . 
is repeated: substituting the detailed asymptotics of G and Go to ( 12. ip and 
comparing the terms in the left and right hand sides, we can arrive at the 
equations 

® i (x) = <$>° i (x)- [ f'^ <5>i(s)ds, xGR 3 , z = 1,...,9, (3.17) 
Jr 3 47r|x — s\ 

the first one being identical to (12. 6p . 

Step 3: completing the proof The estimates |$i(x)| < c(l + |x|) 2 easily 
follow from the integral equations (I3.17P and the evident estimates of the 
same form for Thus, we have $j G Vq. 

As is well known, the polynomials $°j ■ ■ ■ > ^9 fo rm a basis in (unper- 
turbed) Vq. Using the equations (I3.17p . it is not difficult to conclude that 
the same is valid for $ 9 in the perturbed Vq. □ 
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3.3 Proving Theorem [3] 

Recall that we deal with the case m = 2; 

Necessity Assume that a G T 2 holds. It means that there is a function 
h 6 H such that (—A + q) 2 h = in IR 3 \{a} and, in the same time, (—A + q)h 
does not vanish identically. 

Considering (—A + q) 2 h as a distribution on (M 3 ), we see that it is 
supported at x = a. Such a distribution has to be a linear combination of 
the Dirac measure derivatives: 

N N 

(-A + q) 2 h = ^ aj Di5 a = aji-l)® Di6 a =: P N (V a ) 5 a , 
b'l=o lil=o 

where .P/v(V a ) is a differential polynomial of the order N with constant co- 
efficients, which acts on the variable a. Therefore, in accordance with ( 13. 7p . 
the function h has to be of the form 

h{x) = P N {V a ) B{x,a) + h{x) , xeR 3 \{a} 

with a smooth g-biharmonic in M. 3 function h. 
Since 

B(x, a) = c\x — a\ + smoother terms 

near x = a (see ( I3.6|) ). the only option for h to belong to H (to be square- 
summable near x = a) is a>j = for \j\ > 3; hence, we have N < 2, i.e., 

h(x) = P 2 (V a ) B(x, a) + h(x) , x e R 3 \{a}. 

Let us show that h = 0. At first, note that for |x| — > oo the function 
h growths not slower than \x\ 2 . Indeed, otherwise it grows not faster than 
| re | and, by the latter, the function ip := (—A + q ]h is g-harmonic in R 3 and 
tends to zero that is impossible by Convention Hq On the other hand, the 
summand /^(Va) B(x,a) may grow not faster than |x|, whereas the sum h 
belongs to %. Hence, h must vanish identically in M 3 . Thus, we get 

h(x) = P 2 (V ) B(x, a) , xeR 3 \{a}. (3.18) 

12 The possibility tp = is also excluded by the same arguments as were used in the 
proof of Lemma [2] 
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Recall that h G f-L and tend \x\ — > oo. Turning to the structure of the 
asymptotic of B (see (13.81) . ( 13.91) with y = a) and taking into account the 
linear independence of the spherical harmonics (the square brackets in (I3.9P ) , 
we easily conclude that the r.h.s. of (" 13 . 1 8[) can be square-summable if and 
only if the coefficients at the terms, which behave as \x\ and rjr, vanish, i.e., 

P 2 (V ) $<(a) = 0, i = l,... ,9; P 2 (V a ) #(a) = . (3.19) 

The obtained relations for q-harmonic $j are nontrivial only if P 2 (V a ) 7^ 
—A + q(a). Show that it is the case. Indeed, otherwise, we have 

$(a) = (see ( 13~T6D ) = -- P 2 (V)tf(a) = (see (l3~T9D ) = . 

6 

By (12. 3D . we conclude that G( - ,a) £ H, i.e., a is a s-point of the first order, 
i.e., a G T ? that contradicts to the assumption. 

Since $1, ... ,$9 constitute a basis in V q , any g-harmonic polynomial 
p G V q satisfies P 2 (V)p(a) = 0. Hence, [P|] < 8, i.e., the jet degeneration 
at the point a does occur. 

Sufficiency Assume that at a point a G M. 3 the jet degeneration oc- 
curs, i.e., d\ \Vq] < 9. Let C be the 10-dimensional subspace spanned on 
$1, . . . , $9, the degeneration evidently implies d 2 a [C] < 10. By the latter, 
there exists a second order differential polynomial P 2 (V) such that (13. 19ft is 
valid. In the mean time, the same arguments as in the item below (13.191) 
imply P 2 (V) ^ -A + q(a). 

Defining a function h by the r.h.s. of (13.181) and taking into account the 
form of the asymptotic (13.81) and the relations (13.191) . it is easy to see that 
h G U and (-A + qfh = in R 3 \{a}. Hence, a G T q 2 holds. □ 

In the case of a G T q 3 for j > 2, the proof is just more complicated 
in notation. However, it exploits the same idea: a linear combination of 
the relevant g-j-harmonic Green function and its derivatives can belong to 
the space H iff the proper analog of the relations (13.191) does hold that is 
equivalent to a certain jet degeneration at x — a. 

3.4 Radially symmetric potential 

If the potential is of the form q = q(\x\), much more can be said about 
structure of the s-point set T q . As is evident, T q consists of spheres in M 3 
centered at x — 0. Here we briefly announce some results on this case. 
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1. For a fixed I > 0, the set T q l can be characterized as follows. Let ipi be 
a regular solution of the radial Schrodinger equation 

-¥ > f + - L ^ + = °> r>0 (3-20) 

that behaves as <fi(r) ~ r i+1 near r = 0. Note that, in this case, the 
functions 




(\m\ <l + l) 



are g-harmonic and belong to V l q . Introduce the Kram determinants^ 3 ! 



Ap(r):= 

Let iVj be the number of the eigenvalues of the partial Schrodinger 
operator — d 2 r + ^j-^ + g(r) in L 2 (0,oo) and z(A|™) the number of 
zeros of the Kram determinant. The relation 

l— m 

z(AD = ]T(-l) fc AW (3.21) 

fc=0 

is valid. 

2. As we saw in section 2, o"di SC (^) 7^ implies T q ^ 0. Our hypothesis 
is that the converse is also valid, so that the equivalence 

{a disc (H) = 0} & {T q ± 0} (3.22) 

holds. It can be shown that for radial symmetric potentials, if o~^ sc {H) = 
then T q = holds for all k > 1. This fact is established in the fol- 
lowing way. 

One can show that a e IR 3 is a /-jet degeneration point of the q- 
harmonic polynomials iff a certain (I + 1) x (/ + l)-determinant 8i(r) 

It is worthy of noting that this determinants appear in Darboux transform theory 



<fm+i(r) v4+i( r ) ••• V%+i( r ) 
iptir) tp\(r) ■■■ (pf~ m) {r) 
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vanishes as r = \a\. In the mean time, the determinant can be repre- 
sented in the form 

Mr) = A° (AKr)) 2 ...(A|(r)) 2 

that enables one to find the zeros of 81 by the use of (I3.2ip and, hence, 
to verify the presence/absence of the jet degeneration points. 

3. The picture of s-spheres is rather curious. For instance, in the case 
of H a = —A + aq with an interaction constant a > 0, the following 
scenario is possible as a growths: 

(a) for small enough a G [0,«i], one has adi SC (H a ) = and T q k = ill 
for all k > 1 

(b) for a G [cci,^), &disc(H a ) consists of a single eigenvalue A" < 0, 
the eigenfunction Xi being radially symmetric, whereas all T q , k > 1 
are nonempty and consist of the spheres Sk of the radius growing to 
infinity as k — > oo 

(c) Assume that as a passes through a 2 , the second negative eigenvalue 
A2 > A° appears, the eigenfunction x% depending on angle variable as 
Yx(-). In such a case, as a t «2, the certain spheres Sk, k > 2 are 
blown up: lim — 00. 

In the mean time, if A2 is such that X2 * s a l so radially symmetric, 
the picture is quite different: all T*, k > 1 are nonempty and evolve 
regularly as a growths. 

3.5 Comments 

• Some results of the book [TJ on a system of the Schrodinger equations 
f l3.20p with different I > can be interpreted in 's-point terms'. In 
particular, the existence of the zero energy solution of the system, which 
is bounded at r = 0, leads to a certain 'non-standard' factorization of 
the S- matrix that is quite analogous to the effects discussed in sec.l. 

• Let x = a be a s-point of the potential q, and w[ a reversing wave, 
which satisfies 

w{\\ x - a \<\t\ = 0. (3.23) 



24 



The Fourier transform w a ( ■ , k) = (2ir) 2 J R e tkt w a (- ,t) dt obeys 

(-A + q)w a = k 2 w a , 

so that w a ( ■ , k) is a continuous spectrum eigenfunction of the Schrodinger 
operator H. A remarkable fact is that, by ( I3.23p . such an eigenfunc- 
tion is an entire function of the spectral parameter k. In due time, the 
existence of such eigenfunctions was a question, which has been affir- 
matively answered by R.Newton. Now we see that there exists a rich 
family of entire continuous spectrum eigenfunctions {w a ( • , k)}\ a er q as- 
sociated with reversing finite energy waves and parametrized by points 
of the surfaces, which T q consists of. The concrete examples demon- 
strate that these surfaces may be of rather complicated shape (ovaloids, 
toruses, etc). The meaning and role of such eigenfunctions in the scat- 
tering theory are not quite clear yet. 

• The important question whether the set of the finite order s-points 
U^ =1 T g m exhausts T q is still open even in the case of radially symmet- 
ric potentials. Also, the conjecture (I3.22p is not justified yet. 
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